AIRE relies on Z-DNA to flag gene targets for thymic T cell tolerization.

AIRE is an unconventional transcription factor that enhances the expression of thousands of genes in medullary thymic epithelial cells and promotes clonal deletion or phenotypic diversion of self-reactive T cells1-4. The biological logic of AIRE's target specificity remains largely unclear as, in contrast to many transcription factors, it does not bind to a particular DNA sequence motif. Here we implemented two orthogonal approaches to investigate AIRE's cis-regulatory mechanisms: construction of a convolutional neural network and leveraging natural genetic variation through analysis of F1 hybrid mice5. Both approaches nominated Z-DNA and NFE2-MAF as putative positive influences on AIRE's target choices. Genome-wide mapping studies revealed that Z-DNA-forming and NFE2L2-binding motifs were positively associated with the inherent ability of a gene's promoter to generate DNA double-stranded breaks, and promoters showing strong double-stranded break generation were more likely to enter a poised state with accessible chromatin and already-assembled transcriptional machinery. Consequently, AIRE preferentially targets genes with poised promoters. We propose a model in which Z-DNA anchors the AIRE-mediated transcriptional program by enhancing double-stranded break generation and promoter poising. Beyond resolving a long-standing mechanistic conundrum, these findings suggest routes for manipulating T cell tolerance.

Harnessing non-Watson-Crick's base pairing to enhance CRISPR effectors cleavage activities and enable gene editing in mammalian cells.

Genomic DNA of the cyanophage S-2L virus is composed of 2-aminoadenine (Z), thymine (T), guanine (G), and cytosine (C), forming the genetic alphabet ZTGC, which violates Watson-Crick base pairing rules. The Z-base has an extra amino group on the two position that allows the formation of a third hydrogen bond with thymine in DNA strands. Here, we explored and expanded applications of this non-Watson-Crick base pairing in protein expression and gene editing. Both ZTGC-DNA (Z-DNA) and ZUGC-RNA (Z-RNA) produced in vitro show detectable compatibility and can be decoded in mammalian cells, including Homo sapiens cells. Z-crRNA can guide CRISPR-effectors SpCas9 and LbCas12a to cleave specific DNA through non-Watson-Crick base pairing and boost cleavage activities compared to A-crRNA. Z-crRNA can also allow for efficient gene and base editing in human cells. Together, our results help pave the way for potential strategies for optimizing DNA or RNA payloads for gene editing therapeutics and give insights to understanding the natural Z-DNA genome.

Flipons and small RNAs accentuate the asymmetries of pervasive transcription by the reset and sequence-specific microcoding of promoter conformation.

The role of alternate DNA conformations such as Z-DNA in the regulation of transcription is currently underappreciated. These structures are encoded by sequences called flipons, many of which are enriched in promoter and enhancer regions. Through a change in their conformation, flipons provide a tunable mechanism to mechanically reset promoters for the next round of transcription. They act as actuators that capture and release energy to ensure that the turnover of the proteins at promoters is optimized to cell state. Likewise, the single-stranded DNA formed as flipons cycle facilitates the docking of RNAs that are able to microcode promoter conformations and canalize the pervasive transcription commonly observed in metazoan genomes. The strand-specific nature of the interaction between RNA and DNA likely accounts for the known asymmetry of epigenetic marks present on the histone tetramers that pair to form nucleosomes. The role of these supercoil-dependent processes in promoter choice and transcriptional interference is reviewed. The evolutionary implications are examined: the resilience and canalization of flipon-dependent gene regulation is contrasted with the rapid adaptation enabled by the spread of flipon repeats throughout the genome. Overall, the current findings underscore the important role of flipons in modulating the readout of genetic information and how little we know about their biology.

Thermodynamic analysis of Zα domain-nucleic acid interactions.

DNA/RNA molecules adopting the left-handed conformation (Z-form) have been attributed with immunogenic properties. However, their biological role and importance have been a topic of debate for many years. The discovery of Z-DNA/RNA binding domains (Zα domains) in varied proteins that are involved in the innate immune response, such as the interferon inducible form of the RNA editing enzyme ADAR1 (p150), Z-DNA binding protein 1 (ZBP1), the fish kinase PKZ and the poxvirus inhibitor of interferon response E3L, indicates important roles of Z-DNA/RNA in immunity and self/non-self-discrimination. Such Zα domain-containing proteins recognize left-handed Z-DNA/RNA in a conformation-specific manner. Recent studies have implicated these domains in virus recognition. Given these important emerging roles for the Zα domains, it is pivotal to understand the mechanism of recognition of the Z-DNA/Z-RNA by these domains. To this end, we assessed the binding thermodynamics of Zα domain from ORF112 and ADAR1 on T(CG)3 and T(CG)6 oligonucleotides which have high propensity to adopt the Z-conformation. Our study highlights important differences in the mode of oligonucleotide binding by the two Zα domains originating from different proteins. Site-directed mutagenesis was employed together with isothermal titration calorimetry to tease apart finer details of the binding thermodynamics. Our work advances the understanding on binding thermodynamics of Zα domains to their cognate nucleic acid substrates and paves the ground for future efforts to gain a complete appreciation of this process.

The Simple Biology of Flipons and Condensates Enhances the Evolution of Complexity.

The classical genetic code maps nucleotide triplets to amino acids. The associated sequence composition is complex, representing many elaborations during evolution of form and function. Other genomic elements code for the expression and processing of RNA transcripts. However, over 50% of the human genome consists of widely dispersed repetitive sequences. Among these are simple sequence repeats (SSRs), representing a class of flipons, that under physiological conditions, form alternative nucleic acid conformations such as Z-DNA, G4 quartets, I-motifs, and triplexes. Proteins that bind in a structure-specific manner enable the seeding of condensates with the potential to regulate a wide range of biological processes. SSRs also encode the low complexity peptide repeats to patch condensates together, increasing the number of combinations possible. In situations where SSRs are transcribed, SSR-specific, single-stranded binding proteins may further impact condensate formation. Jointly, flipons and patches speed evolution by enhancing the functionality of condensates. Here, the focus is on the selection of SSR flipons and peptide patches that solve for survival under a wide range of environmental contexts, generating complexity with simple parts.

Special Issue: A, B and Z: The Structure, Function and Genetics of Z-DNA and Z-RNA.

It is now difficult to believe that a biological function for the left-handed Z-DNA and Z-RNA conformations was once controversial. The papers in this Special Issue, "Z-DNA and Z-RNA: from Physical Structure to Biological Function", are based on presentations at the ABZ2021 meeting that was held virtually on 19 May 2021 and provide evidence for several biological functions of these structures. The first of its kind, this international conference gathered over 200 scientists from many disciplines to specifically address progress in research involving Z-DNA and Z-RNA. These high-energy left-handed conformers of B-DNA and A-RNA are associated with biological functions and disease outcomes, as evidenced from both mouse and human genetic studies. These alternative structures, referred to as "flipons", form under physiological conditions, regulate type I interferon responses and induce necroptosis during viral infection. They can also stimulate genetic instability, resulting in adaptive evolution and diseases such as cancer. The meeting featured cutting-edge science that was, for the most part, unpublished. We plan for the ABZ meeting to reconvene in 2022.

A widespread pathway for substitution of adenine by diaminopurine in phage genomes.

DNA modifications vary in form and function but generally do not alter Watson-Crick base pairing. Diaminopurine (Z) is an exception because it completely replaces adenine and forms three hydrogen bonds with thymine in cyanophage S-2L genomic DNA. However, the biosynthesis, prevalence, and importance of Z genomes remain unexplored. Here, we report a multienzyme system that supports Z-genome synthesis. We identified dozens of globally widespread phages harboring such enzymes, and we further verified the Z genome in one of these phages, Acinetobacter phage SH-Ab 15497, by using liquid chromatography with ultraviolet and mass spectrometry. The Z genome endows phages with evolutionary advantages for evading the attack of host restriction enzymes, and the characterization of its biosynthetic pathway enables Z-DNA production on a large scale for a diverse range of applications.

Dual conformational recognition by Z-DNA binding protein is important for the B-Z transition process.

Left-handed Z-DNA is radically different from the most common right-handed B-DNA and can be stabilized by interactions with the Zα domain, which is found in a group of proteins, such as human ADAR1 and viral E3L proteins. It is well-known that most Zα domains bind to Z-DNA in a conformation-specific manner and induce rapid B-Z transition in physiological conditions. Although many structural and biochemical studies have identified the detailed interactions between the Zα domain and Z-DNA, little is known about the molecular basis of the B-Z transition process. In this study, we successfully converted the B-Z transition-defective Zα domain, vvZαE3L, into a B-Z converter by improving B-DNA binding ability, suggesting that B-DNA binding is involved in the B-Z transition. In addition, we engineered the canonical B-DNA binding protein GH5 into a Zα-like protein having both Z-DNA binding and B-Z transition activities by introducing Z-DNA interacting residues. Crystal structures of these mutants of vvZαE3L and GH5 complexed with Z-DNA confirmed the significance of conserved Z-DNA binding interactions. Altogether, our results provide molecular insight into how Zα domains obtain unusual conformational specificity and induce the B-Z transition.

ADAR and Immune Silencing in Cancer.

The regulation of immune responses by tumors is central to their survival. By diminishing the production of interferon (IFN) and other inflammatory mediators, tumors enhance immune evasion. Responses initiated by nucleic acid sensors and triggered by dysregulated RNA transcription and cytoplasmic DNA undergo down-modulation in tumors. A protein hub that involves the double-stranded RNA (dsRNA) editing enzyme adenosine deaminase RNA specific (ADAR), the RNase DICER1, and the dsRNA-activated kinase protein activator of PKR (PACT) mediates many of these tumor-intrinsic responses, with in vitro ADAR dependency varying by tumor type (range 11-80%). The central role played by ADAR, both as an enzyme and as a scaffold, sets it as a target for cancer immunotherapy. Therapeutic approaches focusing on the ADAR p150 isoform and its Z-DNA- and Z-RNA-specific Zα domain find support from recent mouse and human studies.

The RNA binding activity of the first identified trypanosome protein with Z-DNA-binding domains.

RNA-binding proteins play a particularly important role in regulating gene expression in trypanosomes. A map of the network of protein complexes in Trypanosoma brucei uncovered an essential protein (Tb927.10.7910) that is postulated to be an RNA-binding protein implicated in the regulation of the mitochondrial post-transcriptional gene regulatory network by its association with proteins that participate in a multi-protein RNA editing complex. However, the mechanism by which this protein interacts with its multiple target transcripts remained unknown. Using sensitive database searches and experimental data, we identify Z-DNA-binding domains in T. brucei in the N- and C-terminal regions of Tb927.10.7910. RNA-binding studies of the wild-type protein, now referred to as RBP7910 (RNA binding protein 7910), and site-directed mutagenesis of residues important for the Z-DNA binding domains show that it preferentially interacts with RNA molecules containing poly(U) and poly(AU)-rich sequences. The interaction of RBP7910 with these regions may be involved in regulation of RNA editing of mitochondrial transcripts.

Oxidative Modification of Guanine in a Potential Z-DNA-Forming Sequence of a Gene Promoter Impacts Gene Expression.

One response to oxidation of guanine (G) to 8-oxo-7,8-dihydroguanine (OG) in a gene promoter is regulation of mRNA expression suggesting an epigenetic-like role for OG. A proposed mechanism involves G oxidation within a potential G-quadruplex-forming sequence (PQS) in the promoter, enabling a structural shift from B-DNA to a G-quadruplex fold (G4). When OG was located in the coding vs template strand, base excision repair led to an on/off transcriptional switch. Herein, a G-rich, potential Z-DNA-forming sequence (PZS) comprised of a d(GC) n repeat was explored to determine whether oxidation in this motif was also a transcriptional switch. Bioinformatic analysis found 1650 PZSs of length >10 nts in the human genome that were overrepresented in promoters and 5'-UTRs. Studies in human cells transfected with a luciferase reporter plasmid in which OG was synthesized in a PZS context in the promoter found that a coding strand OG increased expression and a template strand OG decreased expression. The initial base excision repair product of OG, an abasic site (AP), was also found to yield similar expression changes as OG. Biophysical studies on model Z-DNA strands found OG favored a shift in the equilibrium to Z-DNA from B-DNA, while an AP disrupted Z-DNA to favor a hairpin, placing AP in the loop where it is a poor substrate for the endonuclease APE1. Overall, the impact of OG and AP in a PZS on gene expression was similar to that in a PQS but reduced in magnitude.

Topologically Constrained Formation of Stable Z-DNA from Normal Sequence under Physiological Conditions.

Z-DNA, a left-handed duplex, has been shown to form in vivo and regulate expression of the corresponding gene. However, its biological roles have not been satisfactorily understood, mainly because Z-DNA is easily converted to the thermodynamically favorable B-DNA. Here we present a new idea to form stable Z-DNA under normal physiological conditions and achieve detailed analysis on its fundamental features. Simply by mixing two complementary minicircles of single-stranded DNA with no chemical modification, the hybridization spontaneously induces topological constraint which twines one-half of the double-stranded DNA into stable Z-DNA. The formation of Z-conformation with high stability has been proved by using circular dichroism spectroscopy, Z-DNA-specific antibody binding assay, nuclease digestion, etc. Even at a concentration of MgCl2 as low as 0.5 mM, Z-DNA was successfully obtained, avoiding the use of high salt conditions, limited sequences, ancillary additives, or chemical modifications, criteria which have hampered Z-DNA research. The resultant Z-DNA has the potential to be used as a canonical standard sample in Z-DNA research. By using this approach, further developments of Z-DNA science and its applications become highly promising.

Z-DNA and Z-RNA in human disease.

Left-handed Z-DNA/Z-RNA is bound with high affinity by the Zα domain protein family that includes ADAR (a double-stranded RNA editing enzyme), ZBP1 and viral orthologs regulating innate immunity. Loss-of-function mutations in ADAR p150 allow persistent activation of the interferon system by Alu dsRNAs and are causal for Aicardi-Goutières Syndrome. Heterodimers of ADAR and DICER1 regulate the switch from RNA- to protein-centric immunity. Loss of DICER1 function produces age-related macular degeneration, a different type of Alu-mediated disease. The overlap of Z-forming sites with those for the signal recognition particle likely limits invasion of primate genomes by Alu retrotransposons.

Sequence preference and structural heterogeneity of BZ junctions.

BZ junctions, which connect B-DNA to Z-DNA, are necessary for local transformation of B-DNA to Z-DNA in the genome. However, the limited information on the junction-forming sequences and junction structures has led to a lack of understanding of the structural diversity and sequence preferences of BZ junctions. We determined three crystal structures of BZ junctions with diverse sequences followed by spectroscopic validation of DNA conformation. The structural features of the BZ junctions were well conserved regardless of sequences via the continuous base stacking through B-to-Z DNA with A-T base extrusion. However, the sequence-dependent structural heterogeneity of the junctions was also observed in base step parameters that are correlated with steric constraints imposed during Z-DNA formation. Further, circular dichroism and fluorescence-based analysis of BZ junctions revealed that a base extrusion was only found at the A-T base pair present next to a stable dinucleotide Z-DNA unit. Our findings suggest that Z-DNA formation in the genome is influenced by the sequence preference for BZ junctions.

Hexon and fiber gene changes in an attenuated fowl adenovirus isolate from Malaysia in embryonated chicken eggs and its infectivity in chickens.

Fowl adenovirus (FAdV) is distributed worldwide and causes economic losses in the poultry industry. The objectives of this study were to determine the hexon and fiber gene changes in an attenuated FAdV isolate from Malaysia in specific pathogen-free chicken embryonated eggs (SPF CEE) and its infectivity in commercial broiler chickens. SPF CEE were inoculated with 0.1 mL FAdV inoculum via the chorioallantoic membrane (CAM) for 20 consecutive passages. The isolate at passage 20 (E20), with a virus titer of 108.7TCID50/mL (TCID50, 50% tissue culture infective dose), was inoculated (0.5 mL) into one-day-old commercial broiler chicks either via oral or intraperitoneal routes. The study demonstrated that 100% embryonic mortality was recorded from E2 to E20 with a delayed pattern at E17 onwards. The lesions were confined to the liver and CAM. Substitutions of amino acids in the L1 loop of hexon at positions 49 and 66, and in the knob of fiber at positions 318 and 322 were recorded in the E20 isolate. The isolate belongs to serotype 8b and is non-pathogenic to broiler chickens, but it is able to induce a FAdV antibody titer. It appears that molecular changes in the L1 loop of hexon and the knob of fiber are markers for FAdV infectivity.

C8-Guanine modifications: effect on Z-DNA formation and its role in cancer.

Base modifications are known to affect the structure and function of DNA. C8-guanine adducts from various carcinogenic compounds have been shown to be potent Z-DNA inducers. Hence, it has been hypothesized that Z-DNA plays a role in cancer and other genetic diseases. In this comprehensive review, Z-DNA and the effect of prevalent C8-guanine adducts on the B-Z transition are addressed. The discoveries of Z-DNA binding proteins including ADAR1, E3L, DLM1, and PKZ have suggested the relevance of Z-DNA in living systems. In addition, increasing evidence on the Z-DNA connection to gene transcription and inhibition reveals potential biological functions of the left-handed DNA. Finally, C8-guanine adducts that promote Z-DNA formation can be used as a tool to explore the Z-DNA function and its role in carcinogenesis.

A B-Z junction induced by an A A mismatch in GAC repeats in the gene for cartilage oligomeric matrix protein promotes binding with the hZαADAR1 protein.

GAC repeat expansion from five to seven in the exonic region of the gene for cartilage oligomeric matrix protein (COMP) leads to pseudoachondroplasia, a skeletal abnormality. However, the molecular mechanism by which GAC expansions in the COMP gene lead to skeletal dysplasias is poorly understood. Here we used molecular dynamics simulations, which indicate that an A … A mismatch in a d(GAC)6·d(GAC)6 duplex induces negative supercoiling, leading to a local B-to-Z DNA transition. This transition facilitates the binding of d(GAC)7·d(GAC)7 with the Zα-binding domain of human adenosine deaminase acting on RNA 1 (ADAR1, hZαADAR1), as confirmed by CD, NMR, and microscale thermophoresis studies. The CD results indicated that hZαADAR1 recognizes the zigzag backbone of d(GAC)7·d(GAC)7 at the B-Z junction and subsequently converts it into Z-DNA via the so-called passive mechanism. Molecular dynamics simulations carried out for the modeled hZαADAR1-d(GAC)6d(GAC)6 complex confirmed the retention of previously reported important interactions between the two molecules. These findings suggest that hZαADAR1 binding with the GAC hairpin stem in COMP can lead to a non-genetic, RNA editing-mediated substitution in COMP that may then play a crucial role in the development of pseudoachondroplasia.

Structural and functional studies of a large winged Z-DNA-binding domain of Danio rerio protein kinase PKZ.

The Z-DNA-binding domain of PKZ from zebrafish (Danio rerio; drZαPKZ ) contains the largest ß-wing among known Z-DNA-binding domains. To elucidate the functional implication of the ß-wing, we solved the crystal structure of apo-drZαPKZ . Structural comparison with its Z-DNA-bound form revealed a large conformational change within the ß-wing during Z-DNA binding. Biochemical studies of protein mutants revealed that two basic residues in the ß-wing are responsible for Z-DNA recognition as well as fast B-Z transition. Therefore, the extra basic residues in the ß-wing of drZαPKZ are necessary for the fast B-Z transition activity.

'Z-DNA like' fragments in RNA: a recurring structural motif with implications for folding, RNA/protein recognition and immune response.

Since the work of Alexander Rich, who solved the first Z-DNA crystal structure, we have known that d(CpG) steps can adopt a particular structure that leads to forming left-handed helices. However, it is still largely unrecognized that other sequences can adopt 'left-handed' conformations in DNA and RNA, in double as well as single stranded contexts. These 'Z-like' steps involve the coexistence of several rare structural features: a C2'-endo puckering, a syn nucleotide and a lone pair-π stacking between a ribose O4' atom and a nucleobase. This particular arrangement induces a conformational stress in the RNA backbone, which limits the occurrence of Z-like steps to ≈0.1% of all dinucleotide steps in the PDB. Here, we report over 600 instances of Z-like steps, which are located within r(UNCG) tetraloops but also in small and large RNAs including riboswitches, ribozymes and ribosomes. Given their complexity, Z-like steps are probably associated with slow folding kinetics and once formed could lock a fold through the formation of unique long-range contacts. Proteins involved in immunologic response also specifically recognize/induce these peculiar folds. Thus, characterizing the conformational features of these motifs could be a key to understanding the immune response at a structural level.

Translocation and deletion breakpoints in cancer genomes are associated with potential non-B DNA-forming sequences.

Gross chromosomal rearrangements (including translocations, deletions, insertions and duplications) are a hallmark of cancer genomes and often create oncogenic fusion genes. An obligate step in the generation of such gross rearrangements is the formation of DNA double-strand breaks (DSBs). Since the genomic distribution of rearrangement breakpoints is non-random, intrinsic cellular factors may predispose certain genomic regions to breakage. Notably, certain DNA sequences with the potential to fold into secondary structures [potential non-B DNA structures (PONDS); e.g. triplexes, quadruplexes, hairpin/cruciforms, Z-DNA and single-stranded looped-out structures with implications in DNA replication and transcription] can stimulate the formation of DNA DSBs. Here, we tested the postulate that these DNA sequences might be found at, or in close proximity to, rearrangement breakpoints. By analyzing the distribution of PONDS-forming sequences within ±500 bases of 19 947 translocation and 46 365 sequence-characterized deletion breakpoints in cancer genomes, we find significant association between PONDS-forming repeats and cancer breakpoints. Specifically, (AT)n, (GAA)n and (GAAA)n constitute the most frequent repeats at translocation breakpoints, whereas A-tracts occur preferentially at deletion breakpoints. Translocation breakpoints near PONDS-forming repeats also recur in different individuals and patient tumor samples. Hence, PONDS-forming sequences represent an intrinsic risk factor for genomic rearrangements in cancer genomes.